Asymptotic properties of quantum Markov chains 
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The asymptotic dynamics of quantum Markov chains generated by the most general physically 
relevant quantum operations is investigated. It is shown that it is confined to an attractor space on 
which the resulting quantum Markov chain is diagonalizable. A construction procedure of a basis of 
this attractor space and its associated dual basis is presented. It applies whenever a strictly positive 
quantum state exists which is contracted or left invariant by the generating quantum operation. 
Moreover, algebraic relations between the attractor space and Kraus operators involved in the 
definition of a quantum Markov chain are derived. This construction is not only expected to offer 
significant computational advantages in cases in which the dimension of the Hilbert space is large 
and the dimension of the attractor space is small but it also sheds new light onto the relation between 
the asymptotic dynamics of quantum Markov chains and fixed points of their generating quantum 
operations. 
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I. INTRODUCTION 

Quantum operations, i.e. completely positive and 
trace non-increasing linear transformations acting within 
a Hilbert space, play a central role in quantum theory. 
They describe the most general dynamics of an open 
quantum system which can be induced by unitary trans- 
formations and von Neumann measurements involving an 
additional initially uncorrelated ancillary quantum sys- 
tem [l|. In particular, in quantum information theory 
quantum operations are an indispensable theoretical tool 
for exploring the ultimate limits of quantum theory 0]. 

Recent quantum technological advances Q have stim- 
ulated significant interest in the dynamics of large open 
quantum systems which arc formed by many either in- 
distinguishable or distinguishable elementary quantum 
systems, such as Bose-Einstein condensates [4[ or qubit- 
based quantum networks These technological devel- 
opments raise interesting theoretical questions concern- 
ing the dynamics of large quantum systems under the 
action of iterated quantum operations, so called quan- 
tum Markov chains [f|. This is due to the fact that 
many discrete models in the area of statistical physics by 
which fundamental physical phenomena, such as the ap- 
proach to thermal equilibrium or the asymptotic dynam- 
ics and decoherence of macroscopic quantum systems, are 
explored involve quantum Markov chains. 

For a general classification of quantum Markov chains 
an investigation of their asymptotic dynamics resulting 
from large numbers of iterations of their generating quan- 
tum operations constitutes a natural starting point. In 
this context the natural questions arise which asymptotic 
dynamics is possible for a quantum Markov chain and 
how it is related to spectral properties of its generating 
quantum operation. It is a main intention of this paper 
to address these questions for arbitrary quantum Markov 
chains. 



Recently, some results addressing these questions have 
already been found for special classes of quantum Markov 
chains. In particular, the asymptotic dynamics of 
Markov chains resulting from iterated random unitary 
quantum operations has been investigated in detail. 
These operations are contracting and hence two ma- 
jor results have been established @, Q- Firstly, it has 
been demonstrated that the asymptotic dynamics of such 
a quantum Markov chain is confined to an attractor 
spaces. This attractor space is spanned by all orthog- 
onal cigenspaces of the generating random unitary trans- 
formations which are associated with eigenvalues of unit 
modulus. Based on an analysis of these particular spec- 
tral properties convenient representations for the asymp- 
totic dynamics of this special class of quantum Markov 
chains can be derived. Secondly, it has been shown how 
an orthogonal basis of such an asymptotic attractor space 
is determined by a set of linear equations involving the 
Kraus operators specifying the generating random uni- 
tary quantum operation. Taking advantage of the fun- 
damental contraction property of this special class of 
unital and trace preserving quantum operations recently 
parts of these investigations have been generalized also 
to quantum Markov chains which are generated either 
by trace preserving and unital 0] or by trace preserving 
and sub- unital quantum operations [10( • Although these 
generalizations demonstrate that these quantum opera- 
tions can be diagonalized on their asymptotic attractor 
spaces they still leave important questions open concern- 
ing, for example, the explicit construction of bases of the 
attractor spaces and of their associated dual bases. 

Despite these recent developments it is still unclear 
to what extent similarly powerful results apply to the 
most general and physically relevant quantum Markov 
chains. This is due to the fact that general quantum 
operations are not contracting so that the ideas under- 
lying the proofs of these previous results do not apply. 
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The main purpose of the paper is to close this gap and 
to generalize these previous results to quantum Markov 
chains which are generated by arbitrary quantum opera- 
tions. We present such a generalization revealing that for 
establishing analogous results the contraction property 
is not essential. The desired structural properties can 
be established already from the fact that general quan- 
tum operations are trace non-increasing and completely 
positive. These properties imply the validity of gener- 
alized Schwartz- inequalities [ll|, EH- Saturating these 
inequalities imposes important structural constraints on 
the asymptotic attractor space. It will be shown in the 
following that a saturation of these generalized Schwartz- 
inequalities is possible whenever a strictly positive opera- 
tor exists which is contracted or left invariant by the gen- 
crating quantum operation of a quantum Markov chain. 
Our results also shed new light onto characteristic prop- 
erties of fixed points of quantum operations thereby gen- 
eralizing recent results on the theory of fixed points of 
quantum operations and of noiseless subsystems in quan- 
tum systems with finite-dimensional Hilbert spaces []~3l — 




This paper is organized as follows. In order to estab- 
lish the basic notions in Sec. [TT] known basic proper- 
ties of completely positive quantum maps are summa- 
rized. In Sec. IIIII general properties of the asymptotic 
dynamics of quantum Markov chains generated by gen- 
eral quantum operations are characterized. We show that 
for large numbers of iterations the dynamics is confined to 
an attractor space which is spanned by, in general, non- 
orthogonal (simple) eigenvectors of the generating quan- 
tum operation. Thus, on this attractor space also these 
most general physically relevant quantum Markov chains 
are diagonalizable. The determination of the asymptotic 
dynamics requires the projection of an arbitrary initially 
prepared quantum state onto this attractor space. For 
this purpose one has to construct a dual basis of this at- 
tractor space. In general, this is a complicated task in 
cases in which the dimension of the underlying Hilbert 
space is large even if the dimension of the attractor space 
is small. The construction of these dual bases is ad- 
dressed in Sec. IIVI It is shown that whenever a strictly 
positive linear operator exists which is contracted or left 
invariant by the generating quantum operation it is pos- 
sible to construct this dual basis in a straightforward way 
from a knowledge of the basis of the attractor space. In 
Sec El the relation between the basis of the attractor 
space and Kraus operators defining the generating quan- 
tum operation of a quantum Markov chain is found. Fi- 
nally, in Sec lVIl it is demonstrated that to some extent 
the theoretical treatment of Sees IIVI and [V] can also be 
applied to cases in which the quantum state which is con- 
tracted or left invariant by a quantum operation is not 
strictly positive. 



II. BASIC PROPERTIES OF QUANTUM 
OPERATIONS 

In this section known basic properties of completely 
positive quantum operations are summarized in order to 
introduce our notation and in order to stress properties 
which are important for our subsequent discussion. 

In the following we consider a TV-dimensional Hilbert 
space Jf? equipped with a scalar product (.,.). Let 
B(J%?) be the associated Hilbert space of all linear opera- 
tors acting on Jf? with the Hilbert-Schmidt scalar prod- 
uct (A,B) HS = Ti{A^B}. The corresponding Hilbert- 
Schmidt norm is given by || A ||:= (A, A)hs for all 
A £ B(J$f). Consequently, the induced norm of a linear 
operator S : B(J^) — > B(J^) acting on the space B(,yf) 
can be defined by 

|| S ||= sup \\S(A) || . (1) 
\\M=i 

This latter norm fulfills the important relation (see e.g. 

El) 

IIS HI ^ || (2) 

characterizing a Banach*-algebra [lj|. Thereby, de- 
notes the adjoint map of S with respect to the Hilbert- 
Schmidt scalar product. If a linear map satisfies the re- 
lation || S || < 1 it is called a contraction. 

An arbitrary completely positive linear map V : 
B(J^) — s- B(J^) admits a decomposition into Kraus op- 
erators {A,-}i=i C BpT) [3, [T3], i-e. 

k 

7>(.) = £^(.)4 (3) 

Its adjoint map is also a completely positive map with 
the Kraus operators {A-}*^, i.e. 

k 

pt(.) = E A K)^- ( 4 ) 

3=1 

In our subsequent discussion we call a completely pos- 
itive map V which is also trace non-increasing, i.e. 
^(J) = Ylj AjAj < I, a quantum operation. In the 

special case a map is called 

a quantum channel or a trace preserving quantum op- 
eration. A quantum operation which leaves maximally 
mixed states undisturbed is called unital and satisfies the 
relation V(I) = J2j Aj'^j = ^ ■ Prominent examples of 
unital quantum channels are random unitary operations 
(or random external fields) @, Q ■ In less restrictive cases 
in which V(I) = • AjA^ < I quantum operations are 
called sub-unital. 

In the subsequent sections the dynamics of a quantum 
system with Hilbert space J%? is discussed which is gov- 
erned by the iterated application of a quantum operation 
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V as described by Eq.©. This dynamics constitutes a 
quantum Markov chain with generator V in analogy to 
the corresponding classical case (compare e.g. with Rcfs. 
@ ED])- Thus, starting with an initial quantum state 
< p(0) G B(J$f), after n iterations this state is trans- 
formed into the quantum state p(n) = T jn (p(Q)). Our 
main purpose is to analyze characteristic features of the 
resulting asymptotic behavior of V n and its relation to 
the spectral properties of its generating quantum opera- 
tion V . In general the commutation relation [P,^] = 
need not hold so that in general a diagonalization of the 
generator V is not possible. Therefore, in such cases the 
determination of the resulting asymptotic dynamics of 
V n is complicated as this n-fold iteration may involve 
high powers of non-trivial Jordan normal forms of the 
generator V . 

Despite these possible complications the following use- 
ful theorem holds which captures some basic spectral 
properties of the special class of quantum channels, i.e. of 
trace preserving completely positive quantum operations 

mm- 

Theorem II. 1 IfV is a trace preserving quantum oper- 
ation of the form f3|) and if a denotes the set of all its 
eigenvalues the following statements are equivalent: 

(1) If\eo- then |A| < 1. 

(2) 1 G a. 

(3) For every quantum state < p G B(J4?) the limit 

n— 1 

P= lim -J2V k (p) (5) 

exists and the quantum state < ~p G B(Ji?) is 
a fixed point and eigenvector of V with eigenvalue 
X = 1. 

In our discussion two concepts for analyzing charac- 
teristic properties of quantum operations V play an im- 
portant role, namely eigenspaecs and their ranges. Sup- 
pose A is an eigenvalue of the map V . Its corresponding 
eigenspace is denoted by 

Ker("P — XI) = {X G B(Jt)\V(X) = XX} , (6) 

and the associated range of the map V — XI is denoted 
by 

Ran(P-AJ) = {X G B(Jf)\3Y G B(Jf), X = T{Y) - XY} 

(7) 

Recently, the asymptotic behavior of quantum Markov 
chains which are generated by random unitary trans- 
formations has been investigated in detail 0, H[. As a 
major result it has been shown that all Jordan blocks 
corresponding to eigenvalues A with |A| = 1 are one- 
dimensional and that the asymptotic dynamics of V n can 
be diagonalized on the associated attractor space 

Attr(P) := Ker(P - XI) (8) 

|A| = 1 



formed by the direct sum of all cigenspaces of V with 
eigenvalues |A| = 1. The original idea of the proof of 
this asymptotic structure as presented in Ref. H relies on 
the fact that random unitary transformations are con- 
tractions. Recently, this result has been generalized to 
Markov chains generated by unital quantum channels 
[TJ| or by sub- unital quantum operations [Toj | . Although 
these investigations address a large class of iterated quan- 
tum operations the extension of these results to quantum 
Markov chains generated by arbitrary quantum opera- 
tions, which are generally non-contracting, is still unclear 
as the previously applied arguments do no longer apply. 
In the following we discuss this issue and generalize pre- 
vious results to quantum Markov chains generated by 
arbitrary quantum operations, i.e. completely positive 
and trace non-increasing quantum maps. 

III. ASYMPTOTIC DYNAMICS OF QUANTUM 
MARKOV CHAINS 

In this section the asymptotic dynamics of quantum 
Markov chains generated by arbitrary quantum opera- 
tions is characterized. In general, this characterization 
is complicated by the fact that in general the generating 
quantum operations are not diagonalizable. 

We show that for quantum operations, i.e. trace 
non-increasing completely positive quantum maps V : 
B(Jif) -S- B(Jf ) (compare with Eq.©) with V f {I) < I, 
in the limit n — > oo the dynamics of the iterated map 
V n is confined to an attractor space Attr(V). This lat- 
ter space is spanned by all, in general non-orthogonal 
(simple), eigenvectors of V with eigenvalues A of unit 
modulus, i.e. |A| = 1. Thus, on this asymptotically 
relevant subspace Attr(V) the iterated quantum oper- 
ation V n can be diagonalized. This main result is de- 
rived in two steps which are captured by two main the- 
orems. The first theorem shows that the eigenvalues A 
of a quantum operation V arc confined to the unit cir- 
cle, i.e. |A| < 1, and that the eigenspaces of eigenvalues 
from the set o~\ := {A | |A| = 1} are one dimensional 
but non-orthogonal, in general. The second theorem ex- 
ploits these properties and presents an explicit expression 
for the asymptotic dynamics of an arbitrary initially pre- 
pared quantum state under the influence of the iterated 
quantum operation V n . 

■Theorem III.l For any quantum operation V : 
B(Ji?) — > B(,J%?) as defined by (0), i.e. for any trace 
non-increasing and completely positive map, the follow- 
ing relations hold: 

(1) Any eigenvalue XofV fulfills the relation \ X |< 1. 

(2) The kernel Ker (V — XT) and the range 
Ran (V — XI) of any eigenvalue X with \ X \ =1 
fulfill the relation 

Ker (V - XI) n Ran (V - XI) = {0} . (9) 
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Proof. According to a general theorem [2CJ, [21| a linear 
map V : B(J%°) — > B(Jf?) is completely positive iff the 
quantum state I®V(\&)($>\) resulting from the extended 
mapl<^V : B(Jf) (8 B(JT) -> B(^) <g) B(Jf ) by acting 

on the entangled pure state |4>) = l/V^V^iLi e 
Jf^Jf? is positive, i.e. T<g>7>(|$)($|) > 0. The positivity 
of this quantum state implies 

1 



let us consider an orthonormal basis {\i)} of the Hilbert 
space Jff so that the linear operators form an 

orthonormal basis of the space B(Jff). This implies the 
relation 



N 



(14) 



®V(\®){®\)\j)\k) | =| ^(l\V(\i)(j\)\k) Y< with B Lj = (i\B\j). As the quantum operation T n is 



TV 

^{mm\)\i)^(k\v(mm 



(10) 



with \i)\l), \j)\k) denoting arbitrary elements of an or- 
thonormal basis in the extended Hilbert space Jf? C§> J$f. 
As a consequence, the relations 



N 



E KW>01)I*> 



< 1 



[11) 



Lk=l 



are fulfilled for all possible values of i,j if, in addition, 
the linear map V is trace non-increasing. 

Assuming that A is an eigenvalue of a completely posi- 
tive and trace non-increasing quantum operation V with 
its corresponding eigenvector X we can conclude 



\\V(X) || 2 HA| 2 || A|| 2 = 

N N 

\\ J2 E \k)(l\(k\V(\i)(j\)\l)(i\X\j) f= 

i,j=l l,k=l 
N N 

E i ^wimnmrn i 2 < 

l,k=l »,3=1 
N N 

E \(i\X\j) | 2 e \(k\n\i)(j\)\l) \ 2 <\\* 



(12) 



l.k=l 



with || A ||= ^/Ei!fc=l I ^l A l fc ) I 2 dcnotin g th e Hilbert- 
Schmidt norm of a linear operator A € B(J^). Thus, 
we finally arrive at the first part (1) of this theorem, i.e. 
|A|<1. 

For the second part (2) of this theorem let us assume 
that this statement is false. Thus, there is an operator 
^ A G B(JT) with A G Ker (7> - A/) n Ran (V - XI). 
This implies V{A) = XA and there is an operator 7^ 
B G B(jf ) such that = AS + A. By induction it 

can be verified that 



1 



V n (B) 



1 



X n B 



(13) 



for all numbers of iterations n > 1. As n increases the 
second term on the right-hand side of (|T3"| becomes ar- 
bitrarily small in comparison with the term on the left 
hand side of ([TU)) for |A| = 1. Due to complete positiv- 
ity of the map V this is also true for the first term on 
the right hand side of (fT3"]) , In order to demonstrate this 



completely positive and trace non-increasing we can con- 
clude from Eq. lfTU|) that 



< 



|P n (K)(il)]|fc)l<i 



(15) 



for all n > and for all 1 < k, I < N. So, in view of (TR 
and (fl5]) also all matrix elements of the first term on the 
right-hand side of (|13|) tend to zero for |A| = 1 as n tends 
to infinity. Thus, we finally arrive at the contradiction 
A = to the initial assumption A 7^ 0. □. 

Theorem III. 2 Asymptotically the iterative dynamics 
of any quantum operation V : B(Jff) B(J$?) is given 
by 



Xoo{n) 



E X n X Kl Tr(X x ^X(0)) (16) 

AGCTi ,2—1 



with 



lim 

n— »oo 



X oo (n)~V n {X{0)) 



0. 



(17) 



Thereby, the eigenvectors X\.i are determined by the re- 
lation P(X\,i) = XX\^ (i = 1, ...,d\) with X G <J\ := 
{A||A| = 1}. Their dual vectors X Xi E B(Jt?) (i = 
1, ...,d\) with respect to the Hilbert- Schmidt scalar prod- 
uct are defined by the property Tr(X X ' 1 ' Xy ^1) = 5\\i8ai 
for all A, A' G <J\ and by the property that for X G o~\ each 
X x,z is orthogonal to all eigenspaces Ker(7' — X' I) with 
|A'|<1. 

Proof. Any linear map V : B(J$?) — > B(J4?) can be 
brought into Jordan normal form by an appropriate basis 
transformation T, i.e. 



Tjr 



(18) 



with the (non singular) linear operator T <E B(J%?) defin- 
ing the basis transformation and with J G B(Jff) denot- 
ing the Jordan normal form of V . Thus, J is given by 
a direct sum of Jordan blocks J k of dimensions d k , i.e. 
J = J2k ®Jk- The matrix representation of any of these 
dk x dfc-dimcnsional Jordan blocks is given by 



{Jk)i 



fX k 1 

x k 1 



V 



\ 


1 

x k ) 



(19) 
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Consequently also the iterated map V can be transformed 
to Jordan normal form, i.e. V n = T~ 1 J n T. It is straight- 
forward to demonstrate [22| that the modulus of the 
(i, j)-th matrix element of a Jordan block {Jk) n is given 
by' 

i w i = i ( _ ^ _ ) <i \ k r d * n d K 

Therefore, in the limit of large numbers of iterations n 
the contributions of all Jordan blocks with eigenvalues 
| Afe |< 1 become vanishingly small so that in this limit 
only the Jordan blocks with eigenvalues | Afc |= 1 con- 
tribute to the iterated map V n '. 

Let us now concentrate on these asymptotically con- 
tributing Jordan blocks Jk with eigenvalues | X k |= 1. 
From theorem IIII.ll we can draw the conclusion that 
all these Jordan blocks are one dimensional so that 
the map V n can be diagonalized within the asymp- 
totic subspacc Attr(T') which is spanned by all eigen- 
vectors X\,i, i = l,...,d\ with possibly c?A-fold degen- 
erate eigenvalues | A |= 1. These eigenvectors X\ t i ful- 
fill the relations ^(Xx^) = XX\^ for i = l,...,d\ and 
A G <7i := {A | |A| = 1} and in general they are not 
orthogonal (with respect to the Hilbert-Schmidt scalar 
product). Consequently, in the limit of large numbers 
of iterations n it is sufficient to expand any initial linear 
operator X(0) G B(Jf) in terms of these eigenvectors, 
i.e. 

dx 

A£<7i i— 1 

with Y(0) denoting the part of X(0) which is contained 
in eigenspaces of eigenvalues | A |< 1. The coefficients 
x\_i are given by 

x x ,i = Tr(X x ^X(Q)) (21) 

with X x ^ G B(Jff ) denoting the dual basis vector of X\^ 
which fulfills the relations Th:(X x > l X\i ^) = SwSai for all 
generalized eigenvectors of a Jordan basis, i.e. |A|, |A'| < 
1 and i = l,...,d\, i' = l,...,dy. As V n (Y(0)) tends to 
zero in the limit n -4 oo we finally obtain the result of 
theorem IJTL21 □ 

Theorem IIII.2I generalizes previous results on the 
asymptotic dynamics of quantum Markov chains which 
apply only to subunital channels and unital quantum op- 
erations or even to more restrictive cases, such as unital 
channels, i.e. V{T) = I, V^{I) < I or random unitary 
operations @, Q • It implies considerable simplifications 
as far as the determination of the asymptotic dynam- 
ics is concerned because only eigenspaces correspond- 
ing to eigenvalues of unit modulus contribute and these 
eigenspaces are associated with trivial one dimensional 
Jordan blocks. However, in general these eigenspaces 
spanning the asymptotic attractor space Attr('P) may 
still be non-orthogonal. This complicates the construc- 
tion of the relevant dual vectors which project onto the 



attractor space because for this purpose typically the 
knowledge of the complete Jordan basis of generalized 
eigenvectors of the quantum operation V is required. 
This fact is summarized in the following corollary. 

Corollary III. 3 Let {X\ : {\ with |A| < 1 and i = 
l,...,d\ be a complete Jordan basis of the quantum op- 
eration V . The corresponding non-singular Hermitian 
matrix 

9\i,\>v = Tr{X{jX x , ti ,) (22) 

with |A|,|A'| < 1, i = l,...,d\, i' = l,...,d\> contains all 
relevant information about the non- orthogonality of this 
Jordan basis. The corresponding dual basis with 
\X\ < 1 and i = 1, ...,d\ is then given by 

d x , 

X x ' 1 = J2 Gr%,A'<'*W. (23) 

|A'|<l,i'=l 

In terms of this Jordan basis and its dual the projection 
operator II onto the attractor space Attr(V) is given by 

dx 

n . = X Xii 1r(X x ^ . ). (24) 

A£(Ti -i— 1 

Thus, despite the simplifications resulting from theo- 
rem IIII. 21 in general the determination of the required 
dual vectors of the asymptotic attractor space Attr(T') 
still constitutes a formidable task in particular in cases in 
which the dimension of the Hilbcrt space is large. Never- 
theless, in the subsequent sections it is demonstrated that 
under additional restrictions on the generating quantum 
operations V of a quantum Markov chain both tasks, 
namely the construction of a basis for the asymptotic 
attractor space Attr('P) and the construction of its asso- 
ciated dual basis, can be simplified considerably. 

IV. CONSTRUCTION OF THE DUAL 
ASYMPTOTIC BASIS 

A major open problem which has not been addressed 
in the previous section is whether there exist convenient 
methods which simplify the construction of the dual vec- 
tors X X ' 1 (i = l,...,d\) for all possible eigenvalues 
A G <J\. The knowledge of these dual vectors is crucial 
for projecting any linear operator X(0) or any initially 
prepared quantum state p(0) onto the attractor space 
Attr(-P) according to Ea. (fT^)) . 

In this section it is shown that under the additional 
assumption, that the generating quantum operation V 
of a quantum Markov chain supports the existence of a 
strictly positive quantum state < p G B(Jif) which is 
contracted or left invariant, i.e. V(p) < p, a straightfor- 
ward construction of these dual basis vectors is possible 
from the knowledge of all eigenvectors with eigenvalues 



6 



A £ (Ji . For this construction we exploit the basic prop- 
erty that the generating quantum operation V of a quan- 
tum Markov chain is trace non-increasing and thus fulfills 
characteristic generalized Schwartz- inequalities [IJ EH ■ 
Saturating these inequalities we arrive at the following 
theorem. 

Theorem IV.l If V : B(Jf) -> B(Jf?) is a quantum 
operation with the additional property that there exists 
a quantum state < p £ B(J4?) fulfilling the inequality 
~P(p) < P then for all kernels of eigenvalues A £ o\ the 
following equivalences hold: 

(1) X £ Ker(7> - XI) & Xp- 1 e Ker(7>t - (1/A)I) 

(2) X £ Ker(7> - XI) p~ x X £ Ker(7>t _ (1/A)7) 

(3) X £ Ker{V - XI) & p^Xp £ Ker(7> - A/). 

./Vofe that A £ cti implies the relation A = 1/A iwi/i A 
denoting the complex conjugate of X. 

Proof. For the proof of the first statement (1) we 
investigate the following linear map 

V{X) =V\Xp-^)p^ (25) 

with its adjoint map 

V\X) =V(Xp*)p-*. (26) 

First of all we demonstrate that both maps are contrac- 
tions. Because is subunital, i.e. V^(I) < I, the 
Schwartzoperator inequality 7?tpf)pt(Xt) < pt(XXt) 
applies [ill. IT3[. Thus, for a density operator < p £ 
B(J^) which fulfills the relation 'P(p) < p we obtain the 
inequality 

|| VpO|| 2 = Tr{(V(X))tV(X)} = 

Tv\v\ p -^X^)V\Xp-^) P ^ < 

Tr {p-^X^Xp-^p} =|| X || 2 . (27) 

Consequently || V || = || ||< 1 and both V and are 
contracting linear maps. 

If X £ Ker(7' — XI) a simple calculation reveals that 

\\V\Xp-?)\\ = \\V(X)p-i \\=\X\\\Xp-i\\(28) 

Furthermore, Schwartz's inequality and the contracting 
properties of V and V* imply the inequalities 

|| V^Xp- 1 ' 2 ) W'^Tr^vHXp-^yV^Xp- 1 / 2 )^ = 

Tr^Xp-^ywHXp- 1 ^ < 

|| vv^Xp- 1 / 2 ) mi xp- 1 / 2 1|< 
llvilllvtiiiKXp- 1 / 2 )!! 2 ^!!^- 1 / 2 )!! 2 . (29) 



For A £ (Ti the very left hand side of these inequalities 
equals the very right hand side so that we can conclude 

VVi(Xp-i)=Vi(V(X)p- 1 )^p = Xp-i (30) 

for A £ o\. Thus, using ViX) = XX we finally arrive at 
the relation 

V\Xp- l )= l -Xp- l =XXp- 1 (31) 

in view of |A| = 1. Analogously, it can be demonstrated 
that 

V^V{X P -^) = V (V\Xp- l )p) p- 1 ' 2 = Xp-? (32) 

for A £ <7i, so that we can conclude V(X) = XX provided 
V^Xp- 1 ) = X^Xp- 1 . 

Statement (2) can be proven by applying the same 
reasoning to the linear map W(X) = p 1 ' 2 V^(p~ 1 / 2 X) 
and to its adjoint map W^{X) = p- 1 / 2 V(p 1/2 X). 

Statement (3) is a simple consequence of statements 
(1) and (2). Assuming that X is an eigenvector of V 
with eigenvalue A £ a\, i.e. V(X) = XX, statement (2) 
implies V' ! (p~ 1 X) = i-p^X and statement (1) implies 
V{p~ x Xp) = Xp~ x Xp. □ 

On the basis of this theorem a new scalar product can 
be defined in the space B(Jff). This scalar product is de- 
termined via any strictly positive operator < p £ B(Jt°) 
by (A, B) p = (A, Bp- 1 )ns- It allows us to define the con- 
cept of p-orthogonality by the requirement that two op- 
erators, say A,Be B(Jtff), are p-orthogonal, i.e. A-L P B, 
iff (A, B) p = 0. Based on this concept the following im- 
portant p-orthogonality relations can be proved. 

Theorem IV.2 Let V : B(Jf ) -> B(Jif) be a quan- 
tum operation and let there be a strictly positive operator 
< p £ B(Ji?) such that V(p) < p, then the following 
statements are fulfilled: 

(1) For any eigenvalue X of V with \X\ = 1 kernel and 
range are orthogonal, i.e. 

Ker(V-XI) ± p Ran(P — XI) (33) 

and 

Ker - AI) n Ran - A/) = {0} . (34) 

(2) For any two different eigenvalues X\ and A2 of V 
with I Ai| = |A 2 | = 1 the associated eigenspaces are 
orthogonal, i.e. 

Ker(P-Ai7) _L P Ker(P-A 2 /). (35) 

Proof. Let us consider X £ B(Jf) with X £ KerfP-AI) 
and A £ (Ji so that we obtain the relations (V(X)y = 
T(Xj) = AAt and V^Xp- 1 ) = jXp- 1 from theorem 
IIV.1I Furthermore, let us consider Y £ B(J$?) with Y £ 
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Rar^T 5 — XI), i.e. there exists a / Z e B(Jf) with |A| < 1. According to theorem MI. 21 this latter sub- 
V(Z) — XZ = Y. This implies the relation space does not contribute to the dynamics of the iter- 

ated quantum operation V n in the limit n — > oo. There- 
(X, Y) p = TriX^Yp- 1 ) = Ty{X^V{Z) P ~ 1 ) - (36) fore, Attr(P) is spanned by the p-orthogonal eigenspaces 

XTiiX^Zp- 1 ) = Tn{[P\X p - x )^Z) - ATr(p- 1 X t ^) = Ker ( v ~ A/ ) with A e CT i> Le - 



for A G <j\. 

In order to prove statement (2) let us consider Xi G 
Ker(P - Ail) and X 2 G Ker(V - X 2 I) with Ai ^ A 2 and 
Ai, A2 G o\. This implies the relation 

(X u X 2 ) p = TriXlXzp- 1 ) = ^TriViX^Xzp- 1 ) = 

Xi 

^Tr(XtvHx 2 p- 1 )) = -^TriXlXzp- 1 )) = 

Al A2A1 

-^ r (X 1 ,X 2 ) p . (37) 
A2A1 

Because of Ai, A2 G o\ the relation A 2 Ai = A 2 /Ai applies 
so that (X 1 ,X 2 )p = for Ai ^ X 2 . □ 

Based on these characteristic properties the dual 
asymptotic basis can be constructed in a simple way 
from the knowledge of all eigenspaces Ker("P — XI) for 
all A G <J\. This central result of this section is summa- 
rized in the following theorem. 

Theorem IV.3 Let V : B(Jf) -> B(Jf) be a quantum 
operation with the additional property that there exists 
a strictly positive quantum state < p G B(Jif) with 
T'(p) < P- Under these conditions the dual vectors X X:l 
of the eigenvectors X\,i with X G o~\ and i = l,...,d\ 
which span the asymptotic attractor space Attr('P) are 
given by 

X x >* = Xx,i p-^Trixl.X^p- 1 )}- 1 . (38) 

Proof. The dimensions of the kernel and the range 
of the eigenspace of an arbitrary eigenvalue A fulfill the 
general relation 

Dim(Ker(V- XI)) + Dim(Ran(V - XI)) > Dim(B(Jtf)). 

According to theorem IIV.2I the linear map V fulfills the 
relation Ker("P — XI) _L p Ran(7' — XT) so that we can 
conclude 

Ker(V — XI)) © Ran(V — XI) = B(JT) (39) 

for all A G ci. In addition, theorem IIV.2I also implies 
that Ker(V - Ail) _L p Ker(V - X 2 I) for all Ai ^ X 2 and 
Ai, A2 G o\ which implies the relation 

Ker(7> - XT) © f] Ran{V - XI) = B(Jf°) 

Ae<ri AS<Ji 

with Haso-i Rzn('P — XI) being orthogonal to the asymp- 
totic attractor space Attr('P) and simultaneously contain- 
ing all contributions of eigenspaces Ker(V — XI) with 



Attr(P) = Ker(V - XI) (40) 
A60-1 

and the p-orthogonal projection onto this attractor space 
is achieved by the dual vectors of Eq. ([38|) and by the 
projection operator 

n(.) = J2 X x ,iTr(X x ^ . ). (41) 

Thereby, it is assumed that in case of degeneracy of 
eigenspaces with A G o\ the corresponding cigenstatcs 
X\s are p-orthogonalized. □ 

V. CONSTRUCTION OF THE ASYMPTOTIC 
ATTRACTOR SPACE 

In this section we address the final question how the 
relevant eigenvectors X\^ i = 1, d\ with A G o\ which 
define the asymptotic attractor space Attr('P) are related 
to the Kraus operators which define the generating quan- 
tum operation of a quantum Markov chain. 

For this purpose we use the following property which 
holds for arbitrary completely positive maps ([3]) even if 
they not trace non-increasing. 

Theorem V.l Let V : B{Jif) B(Jif) be a completely 
positive map defined by and let there be a strictly 
positive operator < p\ G B(Jif) satisfying T'(pi) < p\ 
and a positive operator p 2 > satisfying P^{p 2 ) < p 2 , 
then any X G Ker(V — XI) with X G o\ fulfills the set of 
equations 

pzAjXpl 1 = X P2 X P ^A 3 (42) 
for all j G { 1 , . . . , k}. 

Proof. In order to prove this theorem we investigate the 
linear map 

Vj(X) = X^Xp^Aj^pl - ^AjXp- 1 ' 2 (43) 

with X G Ker(7' — XI) and evaluate the following sum of 
non-negative terms 

k 

< ^Tr{^(A)[V/ J (A)]t}=Tr{Apr 1 Atpt (p2) } + 
j'=i 

lApTrjAVr^pxK 1 ^}- 

ATr {Xp^T(X^)p 2 } - ATr {V(X)p^ X^ p 2 } < 

(l-\X\ 2 )Tr{Xp^X^p 2 }. (44) 
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Thereby, the relations T(X) = XX and V(X^) = XX^ 
together with the assumptions V(pi) < Pi and P' (p 2 ) < 
P2 have been taken into account. For A £ <J\ this leads 
to the conclusion Vj(X) = for X £ Ker(V - XI) and 
for each jg{l,..., k}. □ 

Applying this theorem to a trace non-increasing quan- 
tum operation V and to its adjoint (which is generally 
not trace non-increasing) we obtain the following theo- 
rem. 

Theorem V.2 Let V : B(Jf) -> B(Jt) be a quantum 
operation and let there be a strictly positive operator < 
p £ B(J4?) satisfying V (p) < p, then any X £ Ker(7' — XI) 
with X £ o\ fulfills the set of equations 

AjXp- 1 = XXp- x A v 

A)Xp~ x = (l/X)Xp- 1 A], 

Ajp^X = Xp^XAj, 

A]p~ l X = (l/X)p- 1 XA] (45) 

for all j £ {1, . . . , k}. 

Proof. The first equation is a simple application 
of theorem IV. II with p\ = p and p2 — I- The sec- 
ond equation is a consequence of theorem IIV.1I namely 
V^Xp- 1 ) = (l/X)Xp-\ and of theorem EI] applied to 
with pi = I and p2 = p. The third and fourth equa- 
tions are also consequences of these two theorems. □ 

This theorem states that for each eigenvalue A £ o\ 
the set 

Da, p = {x\A j Xp- 1 = XXp- 1 A j ,A]Xp- 1 = XXp- 1 A], 

A jp - l X = Xp^XA^A^p^X = Xp^XA] 
forj£{l,...,fc}} (46) 

includes the eigenspace Ker(7' — XI), i.e. Ker^T 5 — XI) C 
Da, p. In order to address the question under which con- 
ditions we can achieve equality, i.e. Ker("P — XI) = D\ p , 
we add the following two corollaries of theorem lV. 21 which 
state that either the conditions V^(I) < I and V(p) = p 
or the conditions V^(I) = I and V(p) = p are sufficient 
for this purpose. 

Corollary V.3 Let V : B(Jf) -> B(jf?) be a quantum 
operation defined by and let there be a strictly positive 
operator < p £ B(J^) satisfying V(p) = p, then the 
following relations hold: 

(1) Ker(-P - XI) = D A , P for all X £ <r x . 

(2) If X 1 £ Ker(7> - Ai-T), X 2 £ Ker{V - X 2 I) then 
XiX 2 p- x £ Ker(P-AiA 2 J). 

Proof. For the proof of (1) we have to show that 
Da,p Q Ker(7' — XI). For this purpose let us as- 
sume that X £ Da,p and therefore satisfies the rela- 
tions AjXp- 1 = XXp^Aj for j £ {!,..., k}. Multi- 



plying both sides of these equations by pAj and sum- 
ming over all values of j = 1, k we obtain the result 
T(X) = XXp^Tip) = XX, i.e. X £ Ker("P - XI). 

If Xx £ Ker(-P-Aii) and X 2 £ Ker(V-X 2 I) statement 
(1) ensures that X\ £ Da^p and X2 £ Da 2 ,p- Thus, 
from thcorcm llV.TI we can conclude that pA^p -1 £ Da 2jP . 
Consequently, the following two identities arc fulfilled 

A j X 1 X 2 {p- 1 f = A j X 1 p- 1 pX 2 (p- 1 ) 2 

= X 1 X 1 p- 1 A J pX 2 p- 1 p- 1 

= X 1 X 2 X l X 2 (p- 1 ) 2 A j , 

A ] p- 1 X 1 X 2 p- 1 = X 1 p- 1 X 1 A ] X 2P - 1 

= X 1 X 2 p- 1 X 1 X 2 p- 1 (47) 

and analogous equalities for the adjoint linear operators 
Aj so that we can conclude that X\X2p~ 1 £ Da^.p = 
Ker(7> - AiA 2 7). □ 

Unital quantum operations are examples of quantum 
operations fulfilling the assumptions of this theorem. 

This corollary generalizes previous work on the the- 
ory of fixed points and noiseless subsystems of uni- 
tal channels in Hilbert spaces of finite dimensions 
[l5| . Firstly, these results generalize previously devel- 
oped procedures for evaluating asymptotically relevant 
eigenspaces of quantum Markov chains, which apply to 
the special eigenvalue A = 1 only, to all asymptotically 
relevant eigenvalues A £ o\. Secondly, this result applies 
not only to unital channels. Indeed, according to IV. 31 
only one restriction is required for obtaining not only a 
necessary but also a sufficient condition for the construc- 
tion of Ker(A — I) for all A £ 01, namely V(p) = p for 
some strictly positive < p £ B(J^). 

VI. GENERALIZATIONS TO NOT STRICTLY 
POSITIVE QUANTUM STATES 

In our previous considerations we required the ex- 
istence of a strictly positive quantum state such that 
^(p) < P- However, there are quantum operations which 
do not fulfill this condition. Therefore, the natural ques- 
tion arises whether some of our previous results also ap- 
ply to such situations. 

According to theorem III. II a quantum channel V , i.e. 
a trace preserving completely positive map, is always 
equipped with a state p such that V{p) = p. Although 
this state satisfies the condition V(p) < p it need not be 
strictly positive. Let 

p = Y j a i \i> i ){i> i \ (48) 

i 

be its diagonal form with the orthonormal pure states 
\tpi) and a, > 0. The orthogonal projection P p = 
Si l^ili^il onto the support (or range) of this state p 
satisfies the inequality 

v(p p ) < -?—np) < — \ — HP P ) (49) 
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with a m i n (resp. a max ) denoting the minimal (resp. 
maximal) nonzero eigenvalue of the state p. According to 
Ref.[23[ the projection P p always reduces a quantum op- 
eration V, i.e. T(PpXP p ) = P p T(PpXP p )P p . Hence the 
corresponding reduced map V : B(P p Jf) — > B(P p ,J^ a ) 
constitutes a well defined quantum operation. Indeed, 
complete positivity is apparent and if the original map V 
is trace non-increasing (resp. trace preserving) then also 
its reduction is trace non-increasing (resp. trace preserv- 
ing). Moreover, the fixed state p is strictly positive on 
PpJf? '. Hence, all requirements for the applicability of 
our theorems from sections IIIII and IIVI are met for the 
restriction of a quantum channel V onto the subalgebra 
B{P p Jf). 

This construction is feasible for an arbitrary not nec- 
essarily strictly positive state p satisfying V(p) < p. Fur- 
thermore, there is always a maximal state p satisfying 
this property. "Maximal" means that P„ < P p for an 
arbitrary state a satisfying V(a) < a. In this sense the 
orthonormal projection P p on the support of p defines 
the so-called recurrent subspace P p J^f [23[, the maximal 
Hilbcrt space for which our theory applies. 

VII. SUMMARY AND CONCLUSIONS 

We have investigated the asymptotic dynamics of 
quantum Markov chains generated by general quantum 
operations, i.e. by completely positive and trace non- 
increasing linear maps acting in a Hilbert space. It has 
been shown that their resulting asymptotic dynamics is 
confined to attractor spaces spanned by typically non- 
orthogonal eigenspaecs of the generating quantum opera- 
tions associated with eigenvalues of unit modulus. These 
eigenvalues have trivial Jordan blocks so that asymp- 
totically also these most general and physically relevant 
quantum Markov chains are diagonalizable on their at- 
tractor spaces. Furthermore, provided a strictly positive 
operator can be found, which is contracted or left in- 
variant by the generating quantum operation of such a 



quantum Markov chain, an explicit construction of a ba- 
sis of the attractor space and of its associated dual basis 
has been presented. The basis vectors of the attractor 
space are determined by linear equations which depend 
in a simple way on the Kraus operators defining the gen- 
crating quantum operation. Each of these basis vectors 
is related to its dual one by a linear transformation. This 
linear transformation is defined by the inverse of the pos- 
itive operator which is contracted or left invariant by the 
generating quantum operation. This explicit construc- 
tion of the asymptotic dynamics of an arbitrary iterated 
quantum operation is expected to offer significant ad- 
vantages whenever the dimension of the Hilbert space is 
large and at the same time the dimension of the attractor 
space is small. Thus, the theoretical description of the 
asymptotic dynamics described here may be particularly 
useful for applications in large quantum systems. 

Apart from possible practical advantages our discus- 
sion also explicitly demonstrates the close connection be- 
tween the existence of a quantum state which is con- 
tracted or even left invariant by the generating quantum 
operation of a quantum Markov chain and the result- 
ing asymptotic dynamics. This way it generalizes recent 
results on the theory of fixed points of quantum oper- 
ations and of noiseless subsystems in quantum systems 
with finite-dimensional Hilbert spaces |13 - [l5| . 

Finally, we want to emphasize again that the results 
presented here are restricted to finite dimensional Hilbert 
spaces. This is mainly due to arguments involved in the 
proof of theorem IIII.ll They involve a general theorem 
due to Choi [2(| and Jamiolkowski [2l[ which applies to 
finite dimensional Hilbcrt spaces only and which relates 
the complete positivity of a linear map to the positivity 
of an extended map. However, as most arguments in- 
volved in the proofs of our subsequent theorems also ap- 
ply to infinite dimensional Hilbert spaces it is expected 
that at least parts of the construction methods for at- 
tractor spaces presented here can be extended to infinite 
dimensional Hilbcrt spaces. 
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